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POWER AND CHARGE DISSIPATION FROM 
AN ELECTRODYNAMIC TETHER 


The Plasma Motor-Generator project utilizes the influence of 
the geomagnetic field on a conductive tether attached to a LEO 
spacecraft to provide a reversible conversion of orbital energy into 
electrical energy. The behavior of the current into the ionospheric 
plasma under the influence of the geomagnetic field is of significant 
experimental and theoretical interest. Theoretical calculations are 
reviewed which start from Maxwell’s equations and treat the iono- 
spheric plasma as a linear dielectric medium. These calculations 
show a charge emitting tether moving in a magnetic field will 
generate electromagetic waves in the plasma which carry the charge 
in the direction of the magnetic field. The ratio of the tethers 
speed to the ion cyclotron frequency which is about 25 m for a LEO 
is a characteristic length for the phenomena. Whereas for the 
dimensions of the contact plasma much larger than this value the 
waves are the conventional Alfven waves, when the dimensions are 
comparable or smaller, diffraction effects occur similar to those 
associated with Fresnal diffraction in optics. 

The power required to excite these waves for a given tether 
current is used to estimate the impedance associated with this mode 
of charge dissipation. The result, on the order of an ohm, is 
encouraging. 
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I. 


INTRODUCTION 


According to Faraday's law of induction, a voltage will be 
induced in a wire moving across magnetic field lines. Provided there 
is a stationary return path for current, electrical power can be 
extracted (an IxB force will act to reduce the relative speed of the 
wire with respect to the magnetic field), or with the application of a 
reversed voltage greater than the induced voltage, electrical power 
will be expended and propulsion will result. In low Earth orbit when 
the orbital speed, v c is about 8*10 3 m/s, a 10 kilometer long wire 
would have an induced voltage of slightly more than 2 kV. Neglecting 
losses and assuming good electrical contact with the ionosphere, 20 kw 
of power would be equivalent to a propulsion thrust of about 2.5 N.^ 

The essential elements in the proposed system is the establishment 

of a good electrical contact with the ionosphere. The most convenient 

way to establish this contact is with a plasma generator capable of 

producing plasma in such quantities that it will exceed the ionosphere 

(2-5) 

plasma out to a distance of ten meters or so. Hollow cathodes 
will supply both the plasma need for electrical contact with the iono- 
sphere and electrons needed for the electrical current. Electrical 
power/propulsion systems using plasma contactors are generally called 
Plasma Motor/Generators (PMG). 

In the literature, most authors assume Alfven waves are 
responsible for carrying charge from the ends of electrodynamic tethers 
along magnetic field lines into the ionosphere. The theoretical 
calculations are extremely long and the discrepancies of their results. 
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e.g. impedance of the ionosphere vearying from 1 to 10 5 ohms, require 
that careful study of such calculations be undertaken to identify 
results common to all calculation and those effected by model dependent 
assumptions. Consequently, the next section will consider the general 
solutions of Maxwell's equations under the common assumption that the 
ionosphere can be considered linear dielectric material. Subsequent 
sections will consider the effects of tether configuration and 
frequency bands on the predictions. 


II. THE WAVE EQUATION AND ITS SOLUTION 

Maxwell's equations for the rext frame of the ionosphere are 
V- E - *"f> , 

c V'E - - ET , c - £ +«•$ 

when E and B are the electric and magnetic fields, and p and T are 
the total charge and current densities. 

Using the Fourier transform notation 


( 1 ) 


F«.« J '* di 


( 2 ) 


-# ( k * -uW) 

> 

Fcai - (&t J e ,rht '“ l> F(S.w) 

and combining the latter two Maxwell equations yields the wave equation 

r‘k‘ E - *$■*) * (3 ) 


Treating the ionospheric plasma as medium described by a dielectric 
tensor £ : 




( 4 ) 


l 
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and introducing the current density within the plasma by 

$ - *-e = -VCo-e) (5) 

where ' is the plasma conductivity tensor, the wave 

equation can be written as 

c'k'E -c'(t£)~k x uj2 ^* (6) 


where ;f c is the current density within the conducting tether. 

Choosing a coordinate system such that z is parallel to the geomagnetic 


field and x is the direction of tether motion which is assumed for 


simplicity to be perpendicular the geomagnetic field, the dialectric 
.(ID /r. 


tensor of the form 



(7) 


For frequencies in the region of interest, i.e. much less than the 
plasma frequency, ^pe, is very large and E z can be neglected compared 
to D z in the wave equation. 

Without specifying further the frequency dependence of the dia- 
lectric tensor or of the tether current density, the wave equation can 
be solved for the electric field: 

c - l i (h ‘- ) hK * Ck - ; ( t: )/k; - 

t ; " l c‘/l (h’-Mi Xki- ki ) - j*wVt« J 


where k, 2 = k 2 + k 2 and k . 2 = oj 2 ea/c 2 . 
x y m 

As will be explained later, the g 2 - term in the denominator is 
neglible. Consequently the E-field involves Cauchy poles in the 
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complex k-plane at k 2 = k^ 2 , k z 2 = k^ and kj. - 0. 

The power needed for the establishment and radiation of electro- 
magnetic fields and the associated charge and current densities is 
given by Poynting's theorem 


P = W - U 


( 10 ) 


where 


W = - 


1 


j *E d 3 x 


u= s^j 


(E 2 + B 2 )d 3 x 

W is the rate of conversion of mechanical energy to electromagnetic 
energy and U is the field energy. Since the fields do not change with 
time in the rest frame of the tether, the time derivative of U will be 
zero. Further, with the assumption of a cold, collisionless plasma, 
there is no conversion of electromagnetic energy into mechanical and 
the power dissipation takes the form 

j* *E dx 3 (11a) 

J c 

As a result of the tether having a constant speed, \t Qi and a 
constant current, I, both J and E have Fourier transforms of the form 

f(k.u) = F&) &(»-*.«) 


which allows the power dissipation to be recast as the real part of a 


Fourier integral : 

P 


Es j j t *c £)■£(%) 


(lib) 
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Since the current density in the tether, j , is proportional to 

V 

the tether current I, and the electric field is linear in j c> the 
power is proportional to I 2 . Consequently the impedance, Z, of the 
ionosphere required for the maintenance of these charge carrying, 
electromagnetic waves is given by 

Z ~ " iTHF* ] (i 

The charge density in these waves can be calculated from the 
continuity equation 


(° 


OCX. 1) *-V fCK t) or p(k*) = v- k f(k. to\ \/td 


(13) 


where j ( k , 00 ) is the total current density j = j + j = j - (D-E) 


'c u p J c 4 tt 

From the wave equation, it follows that - ^ k*L) = - and thus 


^ A- £(&.*») (14) 

As a result of equations (8), (12) and (14) the electric field, 
the impedance, and the charge density are determined once the frequency 
dependence of the dialectric tensor and the model used for the current 
density of tether are specified. 


III. MATHEMATICAL ORIGIN OF ALFVEN WAVES 

In the frequency regions were the g 2 - term in the denominator of 
equation (8) can be neglected the E-field is seen to have Cauchy-poles 

at k z = k A 2 ’ k2= k A 2 and k x + k y 2 = 0 
where k^ = oj 2 €j./c 2 = k x 2 (^r) *£j.0< x v c ) 

Only the first two of these poles involve k z and thus are candi- 
dates for the generation of waves moving along the direction of the 
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geomagnetic field. As a result, in performing the inverses Fourier 


transform, the most general k z ~integral will be of the form 

/4h 




(15) 


where A and B are functions of k which are non-singular in k z . 
The k 2 = k. 2 pole is really two poles 


k z-' ±k x <: F )£ - r (16) 

which would result in z-dependencies of the form 
When combined with the k x integration this factor would result in 

cx/ , [ ;k,(xi (£) sfzfl. 


If there is to be radiation into the wake of the tether then two 
conditions must be satisfied. First, the appropriate pole is the + 
one of z>o and - om for z<o. The path of integration along the k z axis 
should be appropriately chosen to retain only the physical pole. Second, 
in order that the waves not be damped, i.e. that they are able to carry 
the charge deep into the ionosphere, requires that £i>o. 

This latter condition will be discussed again in the next section where 
it will restrict the frequencies leading to wave generation to three 
regions or bands. 

If the above two conditions are satisfied, then each k - wave 
will have a wake angle (see figure 1) given by 



where v^ = c/e£, is their speed in the z-direction. The k z poles 
remaining to be considered come from k z 2 = k^ - k^ 2 - k 2 . The condition 
that the waves propagate and carry charge along the z-axis, requires 
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-9 


that k* > o, or k x 2 ( (^r-) 2 fi.- 1)> k^ 2 

Since (v c /c) 2 = 7 • 10” 10 , only frequency regions where Cj_ >10 9 

contribute, a severe restriction. Even then, the contributions from the 

k -k plane will be very restricted, 
x y K 

Consequenty, only the pole at the k z 2 = k ft 2 will contribute to 
waves capable of carrying charge into the ionosphere. It should be 
pointed out that contributions to the integral in equation (12) in 
which k z becomes imaginary are themselves imaginary and do not contri- 
bute to the impedance. 


IV. THE DIELECTRIC TENSOR FOR THE IONOSPHERE 

The resonance frequencies for electrons and ions in the ionosphere 
are due to plasma oscillations and cyclotron motion about the B-field. 
The dielectric tensor has the form given in equation (7) with 


- I ~ (18a) 

when R = 1 - Z fs e *)*/(' I * **«»*>) 

l = 1 - X (•£'*)’/( I -/.O.A.) 

2 y _ o J + ions 

and COfk — , lit c t ^ _ electrons 

The sums are to be taken over electrons and all ion species. For a two 


component plasma, 

( to * - uJl ) 0 ~ (fh) ) 

~~ (u»* &•) j 3 ~~ ) 


where = P- * +&* * U V* * *°f e +■ >|(*V -T2, -uj^ )* + 

U 


(18b) 



are the lower and upper hybrid frequencies. A sketch of as a 
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function of co is shown in figure 2. Clearly the requirement that the 
waves not be damped, i.e. e 1 > o^is satisfied in only three frequency 
regions : 

Region I O 6 uj < (19) 

Region II oJ L < (7* 

Region III , , . , 

According to reference 11, numerical values appropriate to an oxygen- 

electron plasma are Cl: 9 3'®' ^ cj* =• J.s-to £.9/0 Mi 

ojft 3.« to 7 Ht , u>« ** 3.5- /o y #t in the first region 

where v^ Q = c(Sl»/w p ^) is the Alfven velocity, and j . 


V. MODELS FOR THE TETHER CURRENT DENSITY 

The simplest model for the current density of the tether is that 
of Estes^. In this model the vertical tether extends from 
-L/jl <_ y <_L/2. and is infinitely thin: 

j y (X,i) - :r Sex ) Sea [ ) - & (/- ) 7 

where x ,= x-v c t,y± = y±L/,z and the 0 function is given by 

9 <*> = [o y < o . 

There is an infinitely thin horizontal tether at each end of the verti- 
cal tether extending from -L /2 <. x # £ L Jx such that the currents 

A A 

split at x' =0 into equal parts, one in the positive x-direction, and 
another in the negative x-direction: 


where x' = x 1 ± L /Z = x - v t ±1/2. Using equation (2) the corres- 

db X C X 

ponding Fourier transformed currents are 
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(20) 


lx - 2S(»-*-*)[*g£7. r - £]*** 


liyl/Z 


K- iSo*-**) 


Thus the Fourier transform of V • ;f*is given by 

k-r* £$<«-*.*) *±A,Ui 


(21a) 


If, as done in reference (7), the horizontal wires are made into sheets, 
the result is 


= -i Z(c 


kjjz . Jti/- 




i J ^ 


It is not possible to tell from the above expressions for k*^ whether 
the vertical wire is infinitely thin or not. Since there is no diver- 
gence of 3* coming from the vertical wire, its cross section does not 
effect T*X It is of course possible to turn the horizontal sheets 
into parallelepipeds. To describe the charge carried into the iono- 
sphere by the plasma cloud generated by a hollow cathode, it might 
suffice to use a spherical or elliptical surface of appropriate size 
over which the divergence of the current density is constant. 

As emphasized by Estes ^ the factor (. S -.!)...,j <x , ^ x ^ ) in the 

k x L x /Z 

expressions for k*j acts as a frequency cut-off for finite L , i.e. the 

X 

dominate frequencies are those for which 

“ * k x v c < * V c /L x (22) 

( 6 ) 

According to Estes , this cut-off for a minimal tether extent (L ='*m) 

X 
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is just below the lower hybrid frequency for the F-layer of the iono- 
sphere. Consequently, in calculating the impedance where the square of 
this cut-off factor occurs, there will be no significant contributions 
from any but the lowest frequency band for the PMG project. 


VI. RESULTS FOR THE IMPEDANCE AND CHARGE DENSITY 


The impedance, Z, of the ionosphere required for the maintenance 
of the charge carrying electromagnetic waves is given by a k-space 
integral of j *»E. Using the electric field found in equation (8) 
yields 


-r* - Py-r, 4ntu> [ [ * j 

j (*)'£&> - —L(kf-*s)kl 


t(k‘ih£ 




( 23 ) 


k^ 2 gives no significant 


As shown in section III, the pole at k 2 
contribution to the generation of waves and to the impedance. The cut- 
off factor discussed in the last section, limits the frequencies to 
those of region I, i.e. cj <_ and where g (jj--)ex Consequently, 
using the simplest model given by equation (21a) yields 


*• A l l J 

mi 


1/2 


( 24 ) 


where k. = k Y (^-) • (1 - (^-) 2 )'^ and k = ^ . With the appropriate 

A X V C 

contour for the k z - integration, the expression for the impedance 


becomes 

z - 


jfe 

ire 


7 ( * (°* . f *A l 

'* - CO 
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The integration over dk y is easily done to give 

• - . — . \2 


z * 


J (fj (' (25> 


where L = L/d, C x = L x /d with d = v c /fti s 25 meters. 

The difference of this final expression for Z and that given in 
reference (6) is just the contribution for j E which that author 

A A 

neglected in his determination of Z. 

If the limit of L < < d < < L is taken, then 


Z -* 2-$( JnZL + Y- /) 


(26) 


where Y is the Euler constant. In this limit the horizontal extent of 
tether is no longer important and agreement is obtained with reference 
7, for their contribution from the first frequency region. Clearly in 
this limit, the physically important frequency cut-off factor coming 
from the finite horizontal extend of the tether system is turned off. 
Consequently, frequencies are no longer limited to the first band and 
the higher bands could give the significant contributions found by 
reference 7 for a horizontal extend of a few tenths of a millimeter, 
i.e. Zj = .35ft, Z n s io 5 ft and 10 4 ft. Such extents are unphysical 

and one should use the results of this reference best as a mathematical 
check. 


A second limit in which L x > > L > > d is perhaps more physical. 


In this limit 


**y*(-L\ 

C Ln) 


(27) 
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first obtained by reference 10. However this limit could not describe 
the PMG project since a plasma having such a large extend in the 
direction of motion would similarly extend in the vertical direction 
and envelope both ends of the tether and thus short-circuit the system. 

The inclusion of the factor of equation (21b), 

due to the extent in the z-direction, brings in a factor of 


kn / i /z \ 

\ MmU/l J 


2 


into the impedance calculation. Since 


k = k (^Oej^ this factor will further strengthen the frequency 

x c 

cut-off. 

The plasma generated by a hollow cathode will itself be distorted 
by the geomagnetic field. Due to the large value of the conductivity 
in the direction of the geomagnetic field compared to that perpendicular 
to the field, one could expect the effective extend of the plasma cloud 
in the field direction, L z , to be very much larger than that in the 
perpendicular direction, L . It might well be that in the PMG project 
the spatial extend of the plasma cloud along the field lines will be 


the critical parameter. 

Turning to the predictions of the models for 
the Alfven wings, equations (8) and (14) yield 



charge density in 
+ OC^Vc^ (28) 


14-16 



Rewriting the denominator of the second term as 






the pole structure remaining, after the non-contributiong pole at 
k 2 = k^ is discarded, is just that encountered in equation (24). The 

result of the k^ and k^ integration is a contribution similar to that 
resulting from the first term in the equation (28) but reduced by a 
factor of (v c /c) 2 . 


The contribution of the first term gives 

where x” = x/d, z' = z (v c /v^) (l-k 2 )"^ 

Only for L > > d does the upper limit of the integration become 

X 

unimportant and diffraction effects due to the (1-k 2 )^ factor relating 
z' and z result in the normal sharp Alfven wings. Figure 3 shows the 
situation for the case of L = 10 m and L = 800m. 

X X 
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VII. CONCLUSION 


This report has attempted to clearly point out the assumptions 
made and common results found in theoretical models put forth to 
explain the dissipation of charge into the ionosphere from the ends of 
a conductive tether. The most dubious of the assumption is that of 
treating the ionosphere in a "cold plasma" dialectric medium. Once the 
assumption is made the spatial extend of the tether need be modeled in 
a reasonable manner. The dimensions of the region over which the 
charge is released into the ionosphere, is shown to play a crucial 
role in restricting the phenomena to low frequencies. If a physically 
reasonable size is used, calculations show that the models predict an 
impedance on the order of a few ohms. 

The adaption of this model to describe the PMG project is not 
easy since the plasma clouds generated at the ends of the tether will 
themselves be strongly affected by the geomagnetic field. Still it 
seems reasonable to accept the impedance value mentioned above as a 
good first approximation. 



REFERENCES 


1. Hite, 6. E., and J. E. McCoy, "Electrodynamic Tether", accepted for 
publication in the American Journal of Physics (1987). 

McCoy, J. E., "Plasma Motor-Generator System Designs for Power and 
Propulsion Interaction", Conference on Tethers in Space, September 
1986. 

McCoy, J.E. , "Electrodynamic Tethers, I. Power Generation in LEO , II. 
Thrust for Propulsion & Power Storage", Intern, Astronautical Fed- 
eration, 35 Congress, Paper 440, (Oct. 1984). 

2. Krishnan, M., R. G. John, W. F. von Jaskowsky, and K. E. Clark, 
"Physical Processes in Hollow Cathodes", AIAA Journal Vol 15 No. 9, 
pp 1217-1225 (1977). 

3. Siegfried, D. E., and P. J. Wilbur, "Studies on an Experimental 
Quartz Tube Hollow Cathode, Electric Propulsion and its Applications 
to Space Missions", Vol 79 of Progress in Astronautics and Aero- 
nautics, pp 262-277 (1981). 

4. Siegfried, D. E., and P. J. Wilbur, "An Investigation of Mercury 
Hollow Cathode Phenomena" , 13th Intern. Electric Propulsion Conf., 
(April 1978). 

5. Siegfried, D. E., "A Phenomenological Model for Orified Hollow 
Cathodes", NASA CR 168026, (Dec. 1982). 

6. Estes, R. D., "Alf.ven Waves from an Electrodynamic Tethered Satel- 
lite System", accepted for publication in J. of Geophy. Res. (1987). 

7. Barnett, A. and S. Olbert, "Radiation of Plasma Waves by a Conduct- 
ing Body Moving through a Magnetized Plasma", J. of Geophys. Res. 

91, 10, 117 (10.135)(1 986 ) . 

8. Dobrowolny, M. and P. Veltri, "MHD Power Radiation by a Large 
Conductor in Motion Through a Magnetoplasma", IL Nuovo Aminto Vol 
9C No. 1, 27-38 (1986). 

9. Rasmussen, C. E., P. M. Banks, and K. J. Harker, "The Exitation of 
Plasma Waves by a Current Source Moving in a Magnetized Plasma: The 
MHD Approximation", J. of Geophys. Res. 90_ No. Al, 505-515 (1985). 

10 Drell, S. D., H. M. Foley, and M. A. Ruderman, "Drag and Propulsion 
of Large Satellites in the Ionosphere: An Alfven Propulsion Engine 

in Space", J. of Geophys. Res. 70 No. 13. 3131-3145 (1965). 

11. Bauer, S. J., Ph ysics and Chemistry in s Pdce» 6 Physics of 
PI anatary Ionosphere, Springer-Verlag p. 135 ( 1973 ) . 


14-20 



